Friedel oscillations of screening in nanotubes 
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PACS 73 . 63 . Fg - Nanotubes 

Abstract. - In 3D and 2D electronic systems the singular contribution to the static permittivity 
e (Kohn singularity) is a small correction to the regular part of e but it results in the leading 
term in asymptotic behavior of the screened potential (Friedel oscillations). In the present letter 
we show that for nanotubes quite different results are valid: e becomes infinitely large at the 
singular point and the Friedel oscillations do not play the dominant role in the screening at the 
large distances. Moreover, the zero and highest cylindrical harmonics of the effective potential are 
screened by quite different mechanisms. 



I Collective effects in quasi- ID systems are remarkable 
by a number of distinguishing features. Probably the most 
non-trivial one relates to the strictly ID system of interact- 
. ing fermions with the linear single particle dispersion law 
' (the Luttinger model). However, there exist also different 
ID objects not described in terms of the Luttinger liquid 
first of all due to necessity to account for the transversal 
subbands and transitions between them. The correspond- 
^ ing example is a nanotube - a hollow cylinder with 2D 
. electron gas on its surface. As usually, availability of the 
' mobile electrons results in renormalization of the electron 
I - electron interaction, i.e., screening. The specific form 

■ of the screened e-e interaction is determined by the effec- 
. five dimension and by the energy spectrum of the electron 
' gas. As is known in 3D plasma the bare Coulomb interac- 
. tion is replaced (due to screening) by the Yukawa law plus 

■ the Friedel oscillations cos {2pFr)/r^ if the electron gas is 
1 degenerate; here pf is the Fermi momentum, r is the dis- 
tance from a point charge, and h — 1. Similar results in 
the 2D case read: the regular part of the screened poten- 
tial at r ^ oo is Qa\/r^ [1], where Q is the initial point 
charge, ob is the effective Bohr radius, whereas the oscil- 
lating contribution is proportional to cos{2ppr — 7r/4)/r^. 
We see that both in 3D and in 2D problems the Kohn 
singularity of e(fc) at fc = 2pp gives the leading term of 
the screened potential for r oo; here e{k) is the static 
permittivity depending on the momentum k. 

Consider now a nanotube with the semiconductor type 
of the single particle energy spectrum: 



Here p is the electron momentum parallel to the tube axis, 
I is the number of the subband, m is the effective mass and 
a is the nanotube radius. Thus, the problem in question 
differs from the similar one for a planar 2D system by 
only quantization of one of the component of momentum: 
Px — l/cL. The considered system is infinite and uniform in 
^-direction (the nanotube axis) and periodic and uniform 
in the azimuthal direction ip . Hence, the electron Green 
function G depends on only differences z — z' , — ip' . In 
the Fourier representation G is diagonal: G{p,p';l,l') — 
Spp'6wG{p,l). 

In what follows we will apply the standard diagram tech- 
nic to find the e — e interaction renormalized by screening 
and we will use the linear theory of screening. We are 
aware that in strictly ID system with linear single-particle 
dispersion law the Luttinger liquid model is valid. The lin- 
ear screening theory relates to only long-wave length limit. 
This theory breaks down for the momentum of the order 
oi Pf just when the Friedel oscillations become essential 
[2-4].i 

However we actually deal with 2D object and we use 
the parabolic dispersion law for electrons. For such a situ- 
ation the Luttinger model is not applicable and we do not 
have any other instrument except RPA to investigate the 
problem. After that it is not surprising that our results 
for the Friedel contribution to the screened potential differ 
qualitatively from the ones obtained in [2-4] and also [5] . 
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^ There is a misprint in the eq. (8) of ref. [2]. Asymptotic be- 
havior of the screened potential U^fj must be l/\x\ In^ rather 
than l/\x\ In otherwise J' U^f f{x)dx diverges while the Fourier 

component U^ff{q) remains finite at g — > 0. 
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The main difference is dependence of our effective 

(scrccincd) potential not only on z (coordinate along the 
tube axis) but also on - azimuthal separation between 
two electrons on the surface of a hollow cylinder. We dis- 
covered that average potential (zero cylindrical harmonic) 
and all the other harmonics are screened qualitatively dif- 
ferent. 

We guess that in the present situation - there is an ex- 
actly solvable model but for only strictly ID system with 
linear dispersion law and an approximate method for gen- 
eral case - it is worth to find general results within RPA 
and to see what experiments will show. This will be done 
in what follows. 

The Gell-Mann-Brueckner [6] chain of electron loops de- 
termining the Fourier component of the screened e-e in- 
teraction V{k, n) comes to the geometrical progression: 



V{k,n) 



FW(fc,n) 



1-Fy(0)(fc,n)n(a;;fc,n)' 
where V^^^ is the bare Coulomb interaction: 



(2) 



L is the number of the highest occupied subband at zero 
temperatiire. pi = \J1'm{Ep — BP) is the Fermi mo- 
mentum of the l-th subband. Eq. (7) demonstrates 
the qualitative difference of the quasi- ID problem and the 
2D and 3D problems. In 2D and 3D systems the Kohn 
singularity a.t k = 2pp gives small corrections to the regu- 
lar part of the polarization operator: ^ In |^| in 3D [7] and 
in 2D case, where £, = {k — '2pf)/'2pf ^ 1- But in 
ID situation no(fc) becomes infinitely large at A; = 2pi and 
the renormalized interaction vanishes at these points. 

To calculate Vo{z) in the regime \z\ >> a wc write in 
the eq. (6) cos {kz) = Re{ex.p{ikz)), and turn the path of 
integration to the upper imaginary semiaxis. At large \z\ 
we have two contributions to Vq{z): Vq from the point k = 
0, where Kq has logarithmic singularity but IIo remains 
finite, and a number of the Kohn contributions from the 
points k = 2pi resulting in the Friedel oscillations Vq{z). 
For the non-oscillating part Vo{z) we get the expansion in 
inverse powers of the value A = In (2|^;|/a) — C ( C is the 
Euler constant) 
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V^°\k,n)=e' [ [ 

J-oo Jo 



ikz—incp 



dzdip 



^Jz'^ + Aa? sin^ (^/2) 

= 47re2/„(|fc|a)if„(|fc|a), (3) 



Here = e^/x, X is the background dielectric constant, 
In,Kn are the modified Bessel functions of the l-th and 
the 3-rd type. The polarization operator (electron loop) 
has a form 



Il{oj;k,n) = —^ E / "^PT- 



fp — k,l — n fp,l 
I — £p-k,l-n — OJ — i5 

(^=+0), 



(4) 



where fp^i = f{sp^i) are the Fermi occupation numbers. 

We are interesting here in the static screening and 
we put from now on = 0. By transforming V{k,n) 
from eq. (2) back to z-space we obtain the expansion 
of the screened interaction in cylindrical harmonics. As 
y(°)(fc, n) and 11(0; k, n) are even functions of n we come 
to the series 

^('^' ^) = / 7^ exp iikz){V{k, 0) + 

OO 

+2^ V{k,n) cos {nip)). (5) 

n=l 

The slowest decreasing term at |2;| — > oo is given by the 
zero harmonic and reads 



Vo{z) = [ —cos{kz) — 
Jo ^ 1 



47re^/o(fca)i^o(fca) 
47re2no(fc)/o(fca)^o(fca) 



where 



no(fc) = n(0;fc,0) = -5-Vln 



k + 2pi 



k-2pi 



, (6) 
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where kq = Uo{k ^ 0) = to[1/po + 2 Ef=i(lM)] A'; 
ttkq is the sum of the partial density of states in occupied 
subbands. Thus, the Coulomb interaction in nanotubes is 
screened rather weakly, Vq{z) ^ e^/(zln^ (^/a))- 

To find the oscillating part Vo{z) we consider a small 
vicinity of the point k = 2pi and note that the factor 
e^IoKo in the integrand of eq. (6) cancels. Thus, we see yet 
another peculiarity of the ID problem: Friedel oscillations 
of the screened potential do not depend on the charge! 
(see below eqs. (9), (18)). This is due to infinitely large 
magnitude of no(fc) &t k = 2pi. Then one can introduce 
the parameter ^ = {k — 2pi)/2pi and take the integral over 
^ from —00 till 00. By the same shift of the integration 
path as described above we arrive at the series in powers 
of l/ln(4p;^;): 



Vo{z) 



2ir'^pi cos{2piz) 

l=-L 



[1- 



2C 



m zln'^{4:piz) ln{4piz) 



(9) 



We see that the amplitude of the Friedel oscillations in 

nanotubes decreases with increasing distance not slower 
than the regular part Vo{z) does. Moreover, the ratio 
Vo/Vo at any z in the order of magnitude can be estimated 
as l/(pFas), and for "metallic" limit of the dense electron 
gas {pFdB ^ 1) the Friedel oscillations become negligibly 
small. 

Qualitatively different results are obtained for non-zero 
harmonics. The factor In{ka)Kn{ka) for n ^ tends to 
the constant l/2n for A; — > and n(A;, n) at = remains 
finite either. Thus, all non-zero harmonics of the Coulomb 
interaction undergo the screening of the dielectric type. 
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that is, their dependence on the distance z coincides with 

that of the bare Coulomb law. Inclcc'cl, from cq. (3) we 
have for the n-th harmonic of the bare interaction 



y(0)(.) = -Q„_V.(l+2^); 



(10) 



V;(°)(z»a):^ 



r(n+l/2) e 



\/7rn! 



z z 



where Qn is the spherical function of the second type. The 
regular part of the screened interaction stems again from 
the region of small fc. Eq. (4) gives: 



Re(n(fc,n)) 



m 



L 

Vln 



Im(n(A;,n)|fe=o = 0. 



and there are two possibilities for the quantity «;„ = 
n(0,n). If the subband occupation for T = terminates 
at / = ±L and ii? > 4L^, then 
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4P' 
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n„-2/n = 7r( 
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(13) 



After that the terms with I = ±Iq in the sum of cq. (12) 
should be replaced by Iln=2io from eq. (13). Thus, for the 
n-th harmonic of the screened potential at large \z\ we 
obtain the contribution from the point k = : 



v^iz) = I>±lZ^!!(«)2.(i + (14) 

VTrn! z z maBti 

The role of the effective dielectric constant for the n-th 
harmonic is played by the quantity 



(1 



1TK 



maBTi 



(15) 



that tends to 1 with increasing n. The dielectric type 
of screening of harmonics with n > 1 can be under- 
stood in terms of the classical electrostatics (suggested by 
M.V.Entin). Each term of the scries in cq. (5) corresponds 
to the one of the expansion in multipoles: n = 1 gives 
dipole-dipole contribution, n = 2 - quadrupole-quadrupole 
and so on. At very large \z\ the fields of all multipoles on 
the nanotube surface are practically normal to the axis, 
hence, they cause only the displacements of electrons in 
the azimuthal direction (p. That is why the z-dependence 
of the e-e interaction does not change; the system simply 



is polarized in accord with the dielectric mechanism when 
only the bound charges arc available. 

As to the oscillating part Vn{z) it is determined by zeros 
of the argument of logarithm in the eq. (11). First of 
all one can sec that such contributions exist not for all 
values of n. For a fixed concentration of electrons ppa < L 
whereas the singularity in the right-hand-side of eq. (11) 
occurs only if pj^a^ > (n — with |Z| < L. Thus, the 
Friedel oscillations of the n-th harmonic exist for n < 2L 
but for n > 21/ we have only non-oscillating contribution 
Vn{z) given by eq. (14). The calculation of the oscillating 
part Vn{z) for n < 2L is totally similar to the derivation 
of eq. (9). The singular points he are given by the relation 



kc=pi± \ p% 



n — I 



(16) 



(11) V V « 

and the polarization operator at A; —> fee takes the form 



Otherwise, for < AL^ and n even the n-th harmonic 
of the potential couples the degenerate states Iq = n/2 
and lo = —n/2. In this case we calculate the limit of the 
uncertainty in eq. (11) at n = 2Zo, A; — > and find 



n w { )ln\ ^' 



(17) 



Ql 



2kcPi{kcPia^ + nl) 



nl{kc-pi) 



The result reads 



K(.) = -— E 



fcc cos(fcc2;) 



mz ^ \vL\zqi) 



(18) 



Thus, only in the case n < 2L the Friedel oscillations 
in nanotubes give a noticeable contribution at large dis- 
tances that is oscillating part of the n-th harmonic of the 
screened potential decreases slower with increasing z as 
compared with the regular part Vn{z). It is worth to note 
also that the Friedel oscillations look like superposition 
of the monochromatic waves with various periods 211 /kc 
(depending both of n and I) rather than the single wave 
cos(2pirr) as it is in 2D and 3D systems. 

The problem of the screened Coulomb potential in nan- 
otubes was considered in [8]. However, the authors have 
given only numerical results for zeroth harmonic Vq{z) 
and did not discuss the Friedel oscillations. 

To conclude, we have analyzed the e-e interaction in 
nanotubes in the frames of the linear screening theory. 
Qualitative difference between quasi-lD systems and 3D 
and 2D systems is established. The Friedel oscillations in 
nanotubes do not determine the asymptotic behaviour of 
the screened potential at large distances in contrast with 
3D and 2D systems. 
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